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ABSTRACT 



Average gravity magnitude and direction approximations allow a 
closed form solution to the equations of motion. The solution yields 
a time-to-go before cutoff as well as a steering function, consisting 
of a thrust attitude and a thrust attitude turning rate. The principles 
of the guidance scheme, which is adaptive for a set of large disturb- 
ances, are outlined with a constant gravity. A spherical earth model 
extension, which includes staging during the guidance phase, is pre- 
sented. Performance data and a comparison of trajectory shaping are 
included where the comparison is against the theoretical optimum 
classical calculus of variations solution. 
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x, y 
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X, 

y 

L 
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X 


tl 


H 

H 

t 2 

R o 


R i 


Rt 


Ki, 

k 2 

m lli m 

*1, 

F 2 

ihi, 

m 2 

T i 

= m 


Definition 

acceleration components in the space fixed x-y system 

velocity components in the space fixed x-y system 

position coordinates in the space fixed x-y system 

velocity components in the injection coordinate system 

position coordinates in the injection coordinate system 

constant thrust attitude angle reference to the £-axis 
required to meet desired velocity end conditions 

thrust attitude angle reference to the x-y axis in the 
space fixed x-y system 

ins tan taneous t ime 

remaining flight time for stage one and stage two, 
respectively 

radius of the earth 

distance from the center of the earth to the vehicle 

distance from the center of the earth to the injection 
point 

coefficients in the steering function 

instantaneous mass of the first and second stage, 
respectively 

thrust of the first and second stage, respectively 

mass flow rate of the first and second stage, respectively 

complete burn-up time of first stage 
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Definition 

complete burn-up time of second stage 
characteristic velocity, first stage 

characteristic velocity, second stage 

average gravity magnitude between the instantaneous 
point and the cutoff point 

total range angle 
average g' c direction 

remaining second stage burn time computed from the 
characteristic velocity equation 

instantaneous velocity deficiency 
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SUMMARY 


An approximate closed form solution of the equations of motion 
allows the derivation of a path adaptive guidance scheme for vehicle 
flight in a vacuum. The scheme is characterized by a limited number 
of presettings and in-flight computation of the guidance parameters. 
The generation of some transcendental functions is required; however, 
no successive approximation procedures are necessary. The guidance 
outputs are time-to-go before cutoff and the steering function which 
consists of a thrust attitude and a thrust attitude turning rate, Kp. 
The turning rate is used for proper attitude control between passes 
through the guidance computer. Since the booster phase of the tra- 
jectories presented are unguided*, perturbations were allowed to build 
up. It was found that the scheme is adaptive for a large set of first 
stage disturbances. The adaptive nature of the scheme also allows it 
to handle second and third stage performance perturbations with only 
a relatively small loss in injection weights when compared to the 
calculus of variation trajectories. No control over the ground range 
from launch to the injection point was attempted. 

SECTION I. INTRODUCTION 

The proposed guidance scheme is based on a steering program 
derived from a set of simplified differential equations of motion. 

The simplification is justified since the implementation of the scheme 
is basically null seeking. The steering program itself has the cal- 
culus of variations as a background. 


^Position and velocity information is not explicitly used for steering. 
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It is well known from literature [1] that, for a flat earth, the 
optimum thrust attitude is given by 


tan (X) 


a 1 + b y t 
c 1 + d f t ‘ 


Law No . 1 


Imposition of orbital conditions without range 

tan (X) = a" + b M t. 

First order expansion of Law No. 1 or 2 yields 

X = a + bt. 


control gives 
Law No. 2 
the form 

Law No . 3 


A comparison survey of spherical earth trajectories [2] using a 
calculus of variations procedure and trajectories using Law No. 3 
has shown that there is little difference in performance. The guidance 
scheme presented in this report uses Law No. 3 as a steering program 
which is updated after each guidance cycle from the state variables 
that can be made available at that time. 


SECTION II. DESCRIPTION 

The principles of the scheme can best be demonstrated by assuming 
a vacuum flat earth with constant gravity, g, as a model. Later in 
the report extensions will be made to a spherical earth. 

The differential equations of motion relative to a vacuum flat 
earth can be written as 


x = a cos X (1) 

y = a sin X - g, (2) 

where a represents the thrust acceleration and X is the thrust 
attitude angle referenced to the x-axis. The symbol, a, is the force 
of the thrust divided by the instantaneous mass of the vehicle. Let 

mi = instantaneous mass of the vehicle, 

m = mass flow rate, 


F = force of the thrust. 
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then, 


or 


a = F/m ! 


a = 


V 

ex 


T - t 


where 


T = m^/rn 


and 


V 

ex 




SP 


F/m. 


The following integrals are given for future use: 


T 



o 


(3) 

(A) 


(5) 



dt 2 = - V 

ex 



T) In 



T 


T 


/ 


a 


o 


tdt 



T 


( 6 ) 


U a tdt 2 = V ( 




where T is the time-to-go from any arbitrary instant 
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The deficiency in velocity required to achieve the desired velocity 
end conditions may be written as 


V. 

i 




n 1 / 2 


- yi + g T )^ 


(9) 


The time-to-go, T, may be determined from equation (9) and the 
characteristic velocity equation 


A V 


V 

ex 



Hence , 


T = T 1 


-AV./Vex^ 

e 


( 10 ) 


Equations (8), (9) and (10) can be solved for T and X for the 
current state variables &x and y x and the required terminal velocity 

components x^ and y^. At this time, assume that T is found from 
equations (9) and (10) through successive approximation methods. 

Later a method will be introduced that eliminates any need for iter- 
ation processes within the scheme. Since the state variables are 
continually changing, the determination of T and the computation of 
X proceed s tepwise using the up-dated values of x x and yi as they are 
obtained. As AVj_ and T approach zero near the end of the 
powered flight, equation (8) becomes indeterminate; however, it will 
be shown later that this difficulty is not serious and can be resolved 
without undue complication. 

Now to enforce an altitude end condition, the parameters Kx and 
K 2 are introduced into Law No. 3: 

X = X - Kx + K 2 t. 


( 11 ) 



6 


cos X « cos X + Kj sin X - K 2 t sin X 

( 12 ) 

sin X ^ sin X - K x cos X + K 2 t cos X • 


Using these trigonometric expressions in equation (2), the following 
integrals are found: 


y T = Yi - g T + 


sin X - Ki cos X 


< 

i — 1 
0 

L 1 J 


L ex - y j 


'V/ 

+ K 2 cos X 




(13) 


= yi - i §t 2 


+ y x T + V 


ex 


sin X - cos X 


(t - T) 



+ V K 2 cos X 
ex 



(T - T) 



(14) 


Now, equations (8), (9) and (10) enforce the desired velocity end 
conditions. The introduction of the parameters K x and K 2 will perturb 
the velocity end conditions if K x and K 2 are not properly selected. 

For the orbital injection case, the first order disturbances are reduced 
in the velocity component normal to the flight path angle at injection. 
The velycity end condition is preserved by setting the difference of 
y^ and y^ to zero. 


y 


T 





K 


cos XV In 
ex 



cos X 



T 


= 0 , 


( 15 ) 



Equation (14), the altitude end condition, and equation (15), the 
preserved velocity end condition, can be solved simultaneously for 
Kjl and K 2 . The equations have the form 
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- A'iKi + B 1 K 2 = 0 (16) 

- + B^K 2 + C 2 = 0, (17) 


where 


A 




1 “ 



(18) 


B 


i 

l 


cos X 


T In 



t j 


(19) 


A 


1 

2 


V 

ex 


cos X 


T 


(t 


T) In 



( 20 ) 


B 


f 


V 

ex 



T 


r 

I T 

L 


(t - T) 



( 21 ) 


C 


t 

2 


yi - y T 



+ YiT + 


V 

ex 


sin X 




( 22 ) 
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Therefore, 


K x = (23) 

MB!, - A^B'o 


and 


K 2 


a i k i 

B' 


(24) 


Thus, the relation X = X - K x + K 2 t can be computed stepwise as the 
current measurements of the state variables are updated. The specified 
presettings are y^, x T and y^. No enforcement of the terminal range is 
attempted. 

SECTION III. SPHERICAL EARTH WITH GUIDANCE OVER TWO STAGES 

The methods employed in the derivation of the guidance equations 
for a spherical earth are essentially the same as those used with a flat 
earth model. However, the equations of motion must be modified to account 
for constantly changing magnitude and direction of the gravity force. 

Since the gravity force field is conservative and the guidance scheme 
is a null seeking system, the change in gravity magnitude and direction 
is approximated. An average gravity magnitude, g% and an average gravity 
direction, 0*, between the current point on the trajectory and the final 
injection point are updated at the beginning of each pass through the 
guidance scheme. Thus, 


R i 2 = x i 2 + ( R o + Yi ) 2 


(25) 


0 X = tan -1 




gi = go 



(26) 


(27) 
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gi + g T 


K - 0i 


( 28 ) 


(29) 

(30) 


where the subscript "1" denotes the instantaneous values and the sub- 
script "T 11 denotes the terminal values. Figure 2 depicts the coordinate 
system used. 

The T) axis of the guidance equations f coordinate system is vertical 
at injection; consequently, the Tj axis must pass through the injection 
point (Fig. 2). Therefore, the | - r\ (guidance) coordinate system requires 
a previous knowledge of the total range angle. The method employed 

for deriving the range angle, 0^, will be shown later. The ^ - Tj 
coordinate system is formed by rotating the space fixed x - y system 
through the range angle, 0^(Fig. 2). The x - y system is considered 
translated to the center of the earth. 


Thus , 


and 



cos 0 
sin 0 T 


- sin # 

' T 

cos 0 




( cos 0rp 
sin 0 T 


- sin 0 T 
cos 0 T 



(31) 


The "t ime -to -go " , T, is computed during each guidance cycle along 
the trajectory as samples of the state variables are taken. The method 
of computing T without the use of successive approximations will be 
shown later. The instantaneous velocity deficiency is defined as 



~ i l ~ g T sin 


/ 


4 - 


\ " 


i)i 


+ 


T cos 0 


( 32 ) 
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FIGURE 2: FLIGHT GEOMETRY 




1] 


The constant thrust attitude, referred to the £ - r) coordinate 
system, necessary to overcome this velocity deficiency at the time T is 


X^ = tan 


.! r " T 


T) 1 + g T cos 0" 


T s in 0' J 


(33) 


For two stages of guidance, the thrust attitude angle must be defined as 

(34) 


Xp = tan 


-i r\ " \ + 8 (Ti + T2) cos & 


It 


T - - g* (Ti + T 2 ) sin 0* j 


where Ti is the time-to-go during the first guidance stage and T 2 is the 
time-to-go during the second guidance stage. T 1 is the time needed to 
deplete the fuel of the first guidance stage at the present mass flow rate. 

Equation (34) enforces the desired cutoff velocity conditions with- 
out any altitude constraint. To enforce the terminal altitude condition, 
it is necessary to introduce the Ki and K 2 parameters into the steering 
program. Since the r\ axis passes through the injection point then the 
terminal al t i tude ^ cond it ion is r) = r) T , likewise, K L and K 2 must be 
chosen such that f) T - fj = 0. As in the flat earth case the steering 
program is 


" Ki + K 2 t. 


It is also required that X| be continuous across staging; therefore, 
during the second stage of guidance, 


*g = Xg.- K x + Ko (T 1 + t). 


(35) 


The equation of motion in the r\ direction is 

.. ~ Vr >v 

T) = asinX^-g cos 0 . 


( 36 ) 


Using the trigonometric approximations given in equations (12), the 
first integral of equation (36) is 
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n T = Ti - g* (Ti + To) cos 0* + 


~ T / T i \ 

sin X. - Ki cos X „ I 1 V In I — ) 

e t ex i \Ti - t l / 


+ V In : 
exn 


— } + K-> cos X ' V T i In ! ^ _ 

T r ..y > - t j ex ^ - yi-, - T- 




+ V T n m 

ex x ' 


\ / 

o) + V e x = Ts in 


M 


- !•; 


(37) 


where 


V is the first stage of guidance, thrust over m, 
exi 

V ex is the second stage of guidance, thrust over m, 

is the initial mass over m of the first stage of guidance, and 
t 2 is the second stage of guidance mass over ni of the second stage. 
The constant thrust attitude X , first integral of equation (36) is 

n T = Ti - g* (Tl + T 2 ) cos 0* + sin X f 


V in 

exi 




- T i 


+ V 

ex2 


In 




(38) 


The first order perturbations caused by the introduction of K-l and Kp 
must be eliminated; hence, the condition f| — rj = 0 yields 


- Ki cos X.. 


V 

L 


In 




+ V 

exi 



Kp COS X. t 


f 

V ex 2 Tl ln Vj 2 - %) T V ex 1 Vg 1 iU Tl - T x " 


, , T i \ 

+ V Ti In r Ti 


+ V Tp ln p ~ 

exp V - Tp - T ' 


- T- 


Y 


( 39) 


0. 
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The second integral of equation (36) is 


nm = m 


T 2 V 


In 


ex i 



+ V 

ex! 



(^1 


Ti) 



+ V 

exp 



( Tp 


T ) In 



+ 


Kp cos X 


I 



V 

exi. 




i 


r 



+ T] V 

exp 


Tp-(T; 


T ) In ( 




Ti) 



V 

exn 


2 



1 . 2 


- t 2 



(T? 



(40) 


Equations (39) and (40) solved for K L and K yield 

B ,C ri 

Kl = Ar-Bj - A r B a 


(41) 


and 


Kp 


A.iKi 


(42) 


where 


Ai 


V 

exi 



In 


(43) 
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+ 




144) 


A2 


cos X^ 



+ V 

ex l 



('Tl 


Ti) 



+ V 


ex 2 


T2 “ ( t 2 


t 2 ) 



(45) 


B 2 


cos X, 




T 2 


V 

ex L 




+ Ti V 


exp 



C 2 = TU - T1 t + T|i ( T 1 + T 2) 


1 * 

J g (Ti 


T 2 ) 2 cos 0 


+ sin X 


I 



In 



(T X - T x ) 


+ V 

ex 2 


T 2 - (t 2 - Tp) In 




(47 ) 


For explanation purposes, a coast period has not be.en included. The 
coast period does not modify the form of the equations; however, the 
coefficients of equations (39) and (40) would have some additional terms. 

Whenever staging occurs, T x is set to zero; thus, no modification 
of the steering program equations is required. Since staging has been 
previously accounted for, the guidance computation proceeds smoothly. 

The problem is finding a proper method of computing 0 and T 
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SECTION IV. THE COMPUTATION OF 2f AND To 

The total range 0 t may be either preset at launch or it may be 
computed in flight. The preset 0^ works very well except for cases 
where the actual range angle exceeds the presetting. If the preset is 
too small, the trajectory tends to be too steep, causing a loss in per- 
formance. When the range angle exceeds the preset 0^, the sign of 0* 
reverses since the injection point no longer lies on the r\ axis. An 
instability may set up in the steering program if this condition prevails. 
To overcome this difficulty, it is necessary to step 0-j- forward as the 
angle approaches 0 T . If a preset is used, some of the adaptivity of 
the scheme is lost for relatively large disturbances. 

The more general cases are covered by computing the total range 
angle in flight. The approximations used in this particular method con- 
tain small errors during the initial portion of guidance; however, these 
errors reduce significantly as the flight progresses. The range angle 
is computed as follows: 

The distance that would be covered by a horizontal flight over a 
flat earth during the first stage of guidance is 


Xi 


V :i T L 


r 

V j T 1 1 _ ( rr t 

ex! , ^ 


Ti) In 


/ 


Ti 


(48) 


Dividing by the terminal radius gives an approximation for the first 
stage range angle. 


- 


-i- 1 

R 


/ 


V-jTi + V 


(Ti - T i) in 


_vn 


- T> , / I 

1 V J 


(49) 


Using equation (49), the second stage velocity deficiency is computed by 


&V. = V- 4 - V ll 
axi 


•ll 


Ti / 
v 


V T - g T ! sin 2 , 


(50) 


where V T is the preset cutoff velocity. The second stage characteristic 
velocity equation is 


f \ 

AV = V 1 n l 7 ] , 

ex? * / 

' Tp - Try 

where T 2 is the second stage estimated time-to-go. 


( 51 ) 
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Equations (50) and (51) yield 


r r 

To = To <1 - EXP I - 


/ Tl \ 

V, + V In ( ! - 


ex 


Li - V m 

<1 Mi - T J T 


V_ - g T i sin 


'01 

i 


1 


'••• L v ex J 

The estimated second stage range angle is computed by 


(52) 


01 


* = b i 


Vi + V In 

exi 


Tl 


- T n 


* . f 0ll \ 

- g Ti Sin [ ~ — ) 


+ V 


ex-: 


Tp - (t - T- ) In 


{ 


\ Tp - Tp 


1 


(53) 


Then the total range angle is 


0rp - 01 + 011 + 012 5 (54) 

and the average gravity direction is 

0* = \ (011 + 012), (55) 

where 0 L is the instantaneous range angle. 

Since 0n has no gravity losses taken into account, 0n will at 
first be too large. However, this error quickly diminishes as the burn 
time decreases. The error in 0 11 causes 0\t> to be too small. The over- 
all effect is to reduce the total error so that the sum of 0n and 0 iy 
produces a surprisingly good estimate of the total range angle Since 

the velocity of the vehicle is continually increasing, the larger portion 
of the flight time takes place in the lower half of the total range angle. 
For particular missions with long burn times, this approximation tends to 
be inaccurate; hence, it is necessary to use a weighted average for 0* 
and g'" . 


The computation of To, without some method of successive approxima- 
tions, requires some knowledge of the length of the second stage time-to- 
go . Either an initial estimate may be preset or the estimated burn time 
from equation (52) may be used. After the first pass through the guidance 
package, the newly computed T ; o is used for the next pass. After the first 
stage burnout and separation occur, the cycle time of the guidance package 
is subtracted from the old Tp and this value used for each succeeding pass 
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Let BT 2 be the error in T 2 , the time-to-go, and let T 2 be the 
estimated time-to-go; then. 


t 2 = t 2 + st 2 . 


(56) 


The velocity deficiency equation may be written as 


AVi = 


P 


i T - 51 - & ( T 1 + T 2 + 6T 2 ) sin ft 


+ 


% “ Ti 


+ g (Ti + + ST^) cos $ 


12 


f ■ 


Let 


(57) 


^ = l r “ ll - g' (T x + T 2 ) sin 0*, 


Al i = h T - hi + g (Tx + Tj>) cos 


(58) 

(59) 


and 


/\ = g 


Aq cos 0 - A | sin / 


* 2 ,* 2 


•J- 2 


(AV ) = (A | ) + (Aq ) . 


(60) 

(61) 


Then 


2 r .* 


(AV t ) = 


A| - g (BTa) sin 0” 


2 r - 2 

iV* •,»- 


Aq + g (BTo) cos 


0*1 , 


(62) 


or 


(AV X ) 2 = (g'‘) 2 (&T 2 ) 2 + 2A (6T 2 ) + (AV') 2 . 


( 63) 


The characteristic velocity equation is 


AV = V In 


e *i \ T i “ T V exs 


+ V In 


t 2 


t 2 - (T 2 + 5T 2 ) J 


(64) 
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But 


t 2 


In 


Now 


Tp - (T's + STs) 


5T; 


= In 


t 2 


t 2 - T 2 


- In 1 - 


&T- 


r 2 - Tp' 


(65) 


In 1 - 


■5T 


2 1 
2 


-5T ; 


+ ... 


t 2 - T y T 2 - T 2 ‘ x t 2 - T^- 

Hence, equation (64) can be written as 


( 66 ) 


AV = V In 
ex x 


. 1_ 
2 


6T 2 


: 1 ” 

2 


+ V In 
ex? 


ex 2 


-6Tp 


T-, - Tp 


+ ... 


(67) 


Assume that T 2 is a reasonably good estimate of the second stage burn 
time; then the terms 1/2 (-ST 2 /t 2 - T 2 ) 2 and higher order are small and 
may be neglected. 

’Let 


L = V In , 

exi Vi - T 


+ V In 
ex 


Tp - Tb 


( 68 ) 


and 


K = 


ex 2 


Tp - T'p 


(69) 


Then, by substituting equations (68) ajid (69) into equation (67) and 
squaring, 5T 2 can be solved from equations (63) and (67), 


(g*) 2 (ST 2 ) 2 + 2A (5T a ) + (AV*) = K 2 (6T 2 ) 2 + 2LK (5T 2 ) + L 2 . 


(70) 


Let 



( 71 ) 
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b = A - LK, 


(72) 


c = (AV^ ) 2 - L 2 . 


(73) 


Then equation (70) can be written as 


-a (oT 2 ) 2 + 2b (6T 2 ) + c = 0. 


(74) 


Therefore, 

5T 2 = L± /b 2 + ac 

a 


(75) 


Since the coefficient b is negative, the positive radical is chosen in 
equation (75) to produce the smaller BT 2 . 


* 


Using T 2 = T 2 + 5T 2 and equation (75) for 0 , X is computed by 


ic 


'X t 


tan 


-i 


* At] + g (ST 2 ) cos 0 
• * * 

-Ag - g (6T 2 ) sin 0 


(76) 


The indeterminate nature of equation (76), mentioned in the description, 
as T 2 approaches zero can be handled in two ways.^ The generally used 
method is to freeze all the guidance parameters, X*, Ki, K 2 and T 2 as 
T 2 < e, where e is some arbitrarily small time-to-go. The data pre- 
sented at the end of the paper was computed using this method. It has 
been found that the constant e is not critical; in fact, an e of up to 
20 seconds causes negligible dispersions in the desired terminal condi- 
tions . 


The second method is to determine the rate of change of the velocity 
components and use the limit equation 





LIM 

•• 1 / 

0 -At}T 2 + 

g cos 0 \ - 

T 2 - >0 

tan 1 ( 

V -AiT 2 - 

* . J>) 

g s in 0 ' J 


(77) 


as prescribed by L'Hospital's Rule. It has been found that this method 
does not appreciably improve the end conditions over the T 2 < e method. 


The steering equation now derived is referenced to the horizontal 
at injection. It is necessary then to rotate the steering function back 
into the space fixed system; hence, 
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X = X g - Kj. - 0 T + K 2 (t - t 0 ), (78) 

where t 0 is the instant in running when X was last computed, and t is 
the running time. 

SECTION V. IMPLEMENTATION PROCEDURE 

The order in which the guidance equations were derived does not 
determine the order of computation. Before computing the steering func- 
tion the total range and the suceeding burning time of each stage must be 
known. Before computing T 2 , the state variables must be rotated into the 
| - T) system of which one axis passes^ through the injection point. T 2 
must now be computed before finding Once X| and T 2 are determined, 

the coefficients of Ki and K 2 can be evaluated. After Ki and K 2 are 
computed, the steering program must then be referenced to the inertial 
coordinate system. A flow diagram will demonstrate the proper computing 
sequence. Figure 3 depicts the flow of two-stage guidance equations in 
the guidance computer. 

Although two stages of guidance are shown, the equations can be 
transformed into one stage simply by setting Ti = 0 and to some 
arbitrary constant. If some engine-out capability is desired, then at 
the instant of the engine failure, the guidance scheme is reinitialized 
by adjusting either and V ex or t 2 and V e x 2 depending on which stage 
is in operation. This reinitialization is a rapid process and has no 
noticeable effect on the overall performance. Engine-out capability is 
not included in the flow diagram. 

The equations needed to compute the guidance constants are 


V? = 


2 

Xl 


y 2 


(79) 


Rf = xf + yf 


(25) 


where the reference system is considered translated to the center of 
the earth. 



gi = go 



(26) 

(27) 
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( 28 ) 



( 29 ) 

( 49 ) 

( 52 ) 


( 53 ) 

( 54 ) 


( 31 ) 


( 58 ) 


( 59 ) 
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FIGURE 3: FLOW DIAGRAM 
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where (■ and r \ ^ are preset values determined by the desired terminal 
path angle. 


A = g 




|^At] 


cos 0 - 


•* * 
sin p 


* 0 • * 0 ~ 
(AV ) 2 = (Ar) ) 2 + (A| ) 2 


L = V 1,(4), V In (- — — -) 

ex i \ T 1 - T V ex 2 V 2 - Tl/ 

V 

ex s 

K = 

t 2 - T 2 


a = K 2 - (g ) 2 


b = A - LK 


c = (AV“) 2 - L 2 


_ b + n/ b 2 + ac 

OT p — " i r . r r r 


T 2 = Tt, + 5T 2 





tan 1 


“Arj + g (ST 2 ) cos 0 - 
-A'l - g" (ST 2 ) sin pT-l 


Ai 


= V 

ex x 


In 



In 



T 


1 


Bi 


V 

ex x 





+ V 

ex 2 


Ti In 



( 60 ) 

( 61 ) 

( 68 ) 

( 69 ) 

( 71 ) 

( 72 ) 

( 73 ) 

( 75 ) 
( 56 ) 

( 76 ) 

( 43 ) 


+ x 



T 


( 44 ) 
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The equations are listed in the form as they appear during the 
derivation. An inspection of this set will reveal that it is possible 
to reduce the computation by combinations. However, the purpose of this 
report is to demonstrate the principles involved in the guidance scheme 
and not to present the scheme in its final form. 

SECTION VI. NUMERICAL INVESTIGATION 

The Saturn V vehicle, suborbital start of the third stage, was 
used as a model to demonstrate the adaptive nature of the guidance 
system presented. However, no data or characteristics of the Saturn V 
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vehicle itself are presented. The tables are all in terms of Delta 
performance and percentages. A realistic comparison was desired; hence, 
for each case presented a calculus of variations program was run for com- 
parison data. This means that the data is not a comparison to some 
nominal standard trajectory, but a comparison to what the calculus of 
variations would have done if faced with the same situation. The guid- 
ance loop was closed after the Saturn V first stage burnout and coast 
period. At this point dispersions were introduced to generate a family 
of orbital injection trajectories. 


The weakest link in the guidance system is the average gravity 
magnitude and direction approximations. For single stages with rela- 
tively short burn times and/or relatively steep trajectories, the approxi- 
mations are excellent. However, as the trajectories become flatter with 
longer burn times, then the larger total range angles cause the approxi- 
mation to begin to breakdown. This is not serious if it is known that 
the missile characteristics will produce such a trajectory. It is only 
necessary to modify the g" equation (29) and the 0 equation (55) to a 
set of weighted average equations. It is even feasible to consider a time 
varying weighted average; however, this is not done in this report. 


To demonstrate the weighted average effects, two cases are presented: 

Case A, /-■■ ~ 2 ( /^ n / ip^) ) (^) 

8 

Case B, 0 = — (Pn + r ) • (80) 


However, no cases are presented with a weighted average gravity magnitude 
g*. Each case was run on a typical Saturn V orbital injection trajectory 
with seven large perturbations introduced to form the family of trajec- 
tories. The results are presented in Table 1. All cases were run with 
the guidance output freezing at To < a method. For this study, an e of 
ten seconds was used. The terminal conditions were so close to the 
desired end conditions that it was not considered necessary to look at 
any smaller c . The results are presented in Table 2. Since it was 
assumed that the vehicle had a perfect autopilot, the errors presented 
are strictly scheme errors. A study of the data presented will show 
that the scheme controls the terminal conditions very tightly with very 
little loss in performance. 


Changing to a weighted average gravity direction causes a change 
in the trajectory shape. Case A mentioned above is an altitude over- 
shoot trajectory with a subsequent performance loss. Although no attempt 
was made to strictly optimize the weighted average, Case B is more nearly 
an optimum trajectory. Figure 4 presents an altitude versus velocity 
plot of Case A, Case B, and the calculus of variations standard case. 

Some trajectory shaping could be accomplished merely by varying the 
weighted average equations (55) and (29). 
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A numerical study was also conducted using an oblate earth model in 
the differential equations of motion. The guidance scheme itself still 
used spherical earth gravity approximations to determine the steering 
coefficients and the time-to-go. This particular study did. not use the 
weighted average gravity magnitude and direction mentioned earlier in 
this section. Since a seventy degree azimuth was used, the oblate earth 
altitude is no longer 185.2 km for a hundred nautical mile orbit refer- 
enced to Cape Canaveral. The first stage was flown with a time poly- 
nomial tilt program with guidance inserted at the ignition of the second 
stage. Table 3 shows the results of this study. The data was generated 
and prepared by the Boeing Company under contract for the Aeroball is tics 
Division of the George C. Marshall Space Flight Center. 


SECTION VII. CONCLUSIONS 

The guidance system outlined in this report is another approach to 
the pa th, adaptive guidance mechanization problem. Once the relatively 
few presettings have been determined, the scheme is capable of handling 
large booster stage disturbances. The scheme requires some subroutines 
like natural logarithm, arc tangent, sine and cosine, square root, etc. 

The guidance output is an initial thrust attitude, a rate of turn 
of that thrust attitude, and a time left to cutoff. At discrete intervals 
of time the state variables sample is updated and a new guidance output is 
generated. Since a turning rate is generated, the discontinuities that 
normally occur in polynomial steering between steps is minimized. An 
indeterminate function in the steering formula occurs if the guidance 
parameters are evaluated very near cutoff. This difficulty is eliminated 
by either applying L T Hospital's rule of limits or by freezing the steering 
constants at some arbitrary time-to-go. The indeterminancy is not serious 
as long as it is accounted for, since it does not affect the terminal con- 
ditions . 

The weakest link in the system is the average gravity direction and 
magnitude computation. However, if it is known what type of trajectory 
the vehicle is required to fly, then a proper weighting function can be 
found that produces an optimum trajectory/ 
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TABLE 1 

PERFORMANCE LOSS OF THE SCHEME IN PERCENT OF C.O.V. MASS IN ORBIT 


Case 

Number 

1 

2 

3 

4 

5 

6 


Initial State Variables 

Performance Loss 

x i 

x i 

*1 

y± 

Percent from Optimum 

Km 

m/ sec 

Km 

m/ sec 

Case A 

Case B 

151.28 

2556.6 

63.974 

930.53 

0. 11% 

0. 04% 

154.00 

2556.6 

63.974 

930.53 

0. 11% 

0 . 04% 

151.28 

2700.0 

63.974 

930.53 

0. 32% 

0. 09% 

151.28 

2500.0 

63.974 

930.53 

0. 07% 

0 . 05% 

151.28 

2556.6 

65.000 

930.53 

0. 11% 

0 . 04% 

151.28 

2556.6 

62.000 

930.53 

0. 11% 

0. 04% 

151.28 

2556.6 

63.974 

1000.00 

0. 14% 

0. 05% 

151.28 

2556.6 

63.974 

850.00 

0. 11% 

0. 03% 


7 
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TABLE II 


TERMINAL CONDITIONS AT CUTOFF OF THE GUIDANCE SCHEME 

DESIRED CONDITIONS ALT. = 185.200 kilometers 
Velocity = 7794.7 m/sec 
Path Angle = 0.0 degrees from the local 
horizontal 


Case 


Case 


A1 titude 


Number 


Kilometers 


1 

2 

3 

4 

5 

6 

7 

8 


185.19995 

185.19991 
185.19987 

185.19997 

185.19998 

185.19997 

185.19992 

185.19998 


Velocity 


Path Angle 


Case B 


A1 titude 


Velocity 


Path Angle 


m/ sec 


7794.679 

7794.680 
7794.642 
7794.686 
7794.680 
7794.679 
7794.675 
7794.695 


Deg. from 
Horizontal 


+.000909 

+.000890 

+.000327 

+.000878 

+.000797 

+.001130 

+.000191 

+.000603 


Kilometers 


185.20000 
185.19999 

185.20001 
185.19999 

185.19998 

185.19999 
185.19999 
185.19999 


m/ sec 

7794.714 
7794.713 
7794.707 

7794.715 
7794.713 

7794.716 
7794.707 
7794.709 


Deg. from 
Horizontal 


+.000275 
+.000268 
+.000165 
+.000424 
+.000246 
+.002540 
+ .002521 
+.000914 
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TABLE 3 

OBLATE EARTH MODEL PERTURBATION EFFECTS ON 
THE GUIDANCE END CONDITIONS 



— 


Resulting End Conditions 

Per turba tion 

Percent 

Deviation 

Alt. 

Km 

Velocity 

M/sec 

Path Angle 
Degrees 

Range 
N. Miles 

Basel ine 
Guidance 

0 

187.58 

7795.40 

89.987 

1425.0 

S-IC Thrust 

- 3 

187.61 

7795.39 

39.988 

1454.1 

S-IC Thrust 

+ 3 

187.58 

7795.39 

89.987 

1414.6 

S-II Thrust 

- 3 

187.62 

7795.38 

89.988 

1462.4 

S-II Thrust 

+ 3 

187.55 

7795.41 

89.985 

1390.9 

S-IVB Thrust 

- 3 


7792.90 

89.989 

1443.2 

S-IVB Thrust 

+ 3 

187.57 

' 

7797.90 

89.982 

1408.1 

S-IC I sp 

| 

187.59 

7795.40 

89.987 

1436.4 

S-IC I sp 

+ 1 

! 

187.57 

7795.40 

89.986 

1413.6 

S-II Isp 

- 1 


I 

7795.40 

89.987 

1434.6 

S-II I sp 

1 

+ I 

.j 

187.58 

7795.40 

89 . 987 

1415.4 

S-IVB I sp 

1 1 

- i | 

187.58 

7795.57 

89.986 

1424.4 

S-IVB I sp 

-f 1 

187.58 

7795.24 | 

89 . 987 

1425.6 

S-IC Wt. 
Uncertainty 

-.3 

187.58 

7795.40 | 

89.987 

1424.2 

S-IC Wt. 
Uncertainty 

t 

+ .3 

187.59 

i 

7795.40 

89.987 

— , , i 
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TABLE 3 (Cont'd) 


Resulting End Conditions 


Perturbation 

Percent 

Deviation 

Alt. 

Km 

Velocity 

M/Sec 

Path Angle 
Degrees 

Range 
N. Miles 

S-II Wt. 
Uncertainty 

-.3 

187.58 

7795.40 

89.987 

1424. 2 

S-II Wt. 
Uncertainty 

+ .3 

187.59 

7795.40 

89.987 

1425.8 

S-IVB Wt. 
Uncertainty 

-.3 

187.58 

7795.40 

89.987 

1424.9 

S-IVB Wt. 
Uncertainty 

+ .3 

187.58 

7795.40 

89.987 

1425.1 

Air Density 

- 3 

187.58 

7795.40 

89.987 

1424.4 

Air Density 

+ 3 

187.58 

7795.40 

89.987 

1424.4 

Drag 

-10 

187.58 

7795.41 

89.987 

1422.2 

Drag 

+10 

187.59 

7795.40 

89.987 

1428.2 























































32 


REFERENCES 


1. Fried, B. D. , "On the Powered Flight Trajectory of an Earth Satellite," 
Jet Propulsion, Vol. 27, June 1957, pp. 641-643. 

2. Lawden, D. F., "Optimal Rocket Trajectories," Jet Propulsion, Vol. 27, 
December 1957, p. 1263. 

3. Chandler, D. C. , H. Horn, and D. T. Martin, "An Iterative Guidance 
Scheme and Its Application to Lunar Landing," NASA Tech. Note in 
preparation. Presently Marshall Space Flight Center report, 
MTP-AERO-63-1 1 (limited distribution only). 



